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This contribution summarizes some recent work on gravitational-wave astronomy 
and, especially, on the generation and detection of relic gravitational waves. We 
begin with a brief discussion of astrophysical sources of gravitational waves that 
are likely to be detected first by the coming in operation laser interferometers, 
such as GEO, LIGO, VIRGO. Then, we proceed to relic gravitational waves em- 
phasizing their quantum-mechanical origin and the inevitability of their existence. 
Combining the theory with available observations, we discuss the prospects of di- 
rect detection of relic gravitational waves. A considerable part of the paper is 
devoted to comparison of relic gravitational waves with the density perturbations 
of quantum-mechanical origin. It is shown how the phenomenon of squeezing of 
quantum-mechanically generated cosmological perturbations manifests itself in the 
periodic structures of the metric power spectra and in the oscillatory behaviour 
of the CMBR multipoles C; as a function of I. The cosmological importance of 
the theoretically calculated (statistical) dipole moment Ci is stressed. The paper 
contains also some comments on the damage to gravitational-wave research in- 
flicted by the "standard inflationary result" . We conclude with the (now common) 
remarks on the great scientific importance of the continuing effort to observe relic 
gravitational waves, directly or indirectly. 

1 Which sources of gravitational waves will be detected first ? 

The comprehensive contributions of Coccia and Giazotto i have described 
the outstanding experimental effort that is now going on at detecting gravita- 
tional waves. As we all know, great expectations are related with the coming 
in operation laser interferometers: the British-German GEO6OQ-J2, the two 
American instruments LIGO □, and the French-Italian VIRGO It is im- 
portant to know in advance, and to be prepared to monitor, those sources 
of gravitational waves that are likely to be detected first in the ongoing and 
forthcoming observations. Before proceeding to the main topic of my contri- 
bution - relic gravitational waves - I would like to present a theorist's view 
on this issue. I will briefly report on the conclusions of a detailed study un- 
dertaken in Ref.Q. By comparing our theoretical analysis with only 3 out of 
several experimental programs, we, of course, do not mean to diminish the 
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importance of work of other experimental groups. 

It is common to divide possible sources of gravitational waves in three 
broad categories: quasi-periodic, explosive, and stochastic. We are aware of 
definitely existing astrophysical sources belonging to each category. To be 
interesting from the point of view of its detection, the source should be suffi- 
ciently powerful, should fall in the frequency band of the detector, and should 
occur reasonably often during the life-time of the instrument. After having 
analysed many_possible sources, and taking all the factors into account, the 
authors of Ref.tj gave their preference to compact binaries consisting of black 
holes and neutron stars. It is clear that in order to radiate large- intensity 
gravitational waves at frequencies accesible to ground-based instruments, the 
objects forming a pair should be massive and should orbit each other at very 
small separations - a few hundreds kilometers. According to the existing 
views, such massive objects can only be the end-products of stellar evolution 
- neutron stars (NS) and black holes (BH). Compact binary systems of our 
interest are coalescing pairs that are only tens of minutes away from the final 
merger, i.e., from the formation of a resulting black hole or, possibly, from 
another spectacular event, such as a gamma-ray burst. The important ques- 
tion is how many such close binaries exist in our Galaxy and at cosmological 
distances. This determines the event rate, that is, the number of coalescence 
events that are likely to occur in a given cosmological volume during, say, 
a 1 year of observations. The event rate is partially constrained by the ex- 
isting observations of binary pulsars, but its evaluation mostly relies on the 
numerical modelling of diverse evolutionary tracks of massive binary stars. 

The event rates are sensitive to several evolutionary parameters, one of 
which is the so-called kick velocity parameter wq. This parameter enters the 
distribution function for velocities that may be imparted on a neutron star 
or a black hole at the moment of its formation. It is believed that wq lies 
in the interval (200 - A00)km/sec. In Fig. | we plot the NS+NS, NS+BH, 
and BH+BH merging rates as functions of Wq. This calculation has been per- 
formed for a typical spiral galaxy of the mass lO n M0 and under a reasonable 
choice of other evolutionary parameters 0. 

It follows from this graph that the NS+NS rate is expected to be at the 
level 3 x XQ~ 5 year~ l , whereas the BH+BH rate is, at least, one order of 
magnitude lower, i.e., 3 x 10~ e year~ 1 . These rates for a typical galaxy, TZg, 
define the rates for a given cosmological volume, IZv, which includes many 
galaxies. When deriving IZv, we have used a conservative estimate for the 
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Figure 1. NS+NS, BH+NS, and BH+BH merging rates in a 10 11 Mq galaxy as functions 
of the kick velocity parameter wo for the Lyne-Lorimer kick velocity distribution. 



baryon content of the Universe, and arrived at the relationship 



Thus, within the volume of radius r — 100 Mpc, and during 1 year, one 
expects that there will be 3 of the NS+NS coalescences and only 0.3 of the 
NS+BH or BH+BH coalescences. The increase of the radius r to r = 200 Mpc 
increases the volume and the event rates by the factor 8. 

The calculated event rates and gravitational wave amplitudes emitted by 
individual binaries should be compared with the instrumental sensitivities. 
The important fact is that the mass of a typical neutron star is 1.4 M Q , 
whereas the mass of a typical stellar black hole is (10 - 15) M Q . Thus, a 
pair of black holes is a more powerful source of gravitational waves than a 
pair of neutron stars. The black hole binaries are less numerous than the 
neutron star binaries, but they can be seen, by a given instrument, from 
much larger distances. The observationally required situation is when the 
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Figure 2. Signal-to-noise ratio, in initial interferometers, as a function of total mass, for 
inspiral signals from binaries of equal masses at 100 Mpc and averaged over source inclina- 
tion. 



detector's sensitivity allows one to see a few events per year, and at the signal 
to noise ratio S/N level of, at least, 2. Thus, the problem now is to find out 
the S/N for binaries of pa, given total mass M. 

In Fig. U we plot □ the (S/N) as a function of M. We have used the 
sensitivity curves of the initial interferometers in LIGO, VIRGO and GEO. 
It is assumed that the binaries consist of equal mass components, so the total 
mass of a NS+NS system is around 2.8M , while the total mass of a BH+BH 
system is around (20 - 30) Mq. The S/N curves are shown for binaries placed 
at r = 100 Mpc. The increase of the distance to r = 200 Mpc decreases the 
values of the plotted functions by the factor 2. 

It is seen from this graph that the S/N is less then 1 for the NS+NS 
binaries. Although the required 3 mergers of neutron stars do occur in the 
volume with the radius r = 100 Mpc, these sources cannot be regarded as 
detectable by the initial interferometers. On the other hand, for a pair of 
black holes with the total mass M — (20 — 30)Af Q , the S/N can be at the 
level of 4, or even larger (in the LIGO detectors). The number of the BH+BH 
events in the discussed volume is too small, though. However, the S/N is still 
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at the level of 2, or better than 2, if the binary black hole is placed at the twice 
as large distance, i.e. r = 200 Mpc. As was shown above, one can expect 
about 2-3 BH+BH events in this larger volume. This is why we believe U 
that the coalescing black holes will probably be the first sources detected 
by the initial laser interferometers. It appears that other possible sources of 
gravitational waves are not likely to be detected at all by these instruments. 
The situation will become much better, and for all known sources, when the 
advanced interferometers become operational. And, of course, there exists a 
nonzero probability that the first detected signals wilLbe totally unexpected 
to everyone, or something what the authors of Ref. Li did not even include 
in their list of options. For additional information on astrophysical sources of 
gravitational waves and their detection, see the comprehensive contributions 
by Ferrari El and Vallisneri El. 

Relic gravitational waves are not expected to be accessible to the ini- 
tial laser interferometers. The prospects are much better, though, for the 
advanced ground-based instruments and for the space-based interferometer 
LISA El. At very large scales, relic gravitational waves should also manifest 
themselves through the anisotropy and polarisation of the cosmic microwave 
background radiation (CMBR) . It is fundamentally important to discover and 
study relic gravitational waves. They are the only clue to the processes that 
have taken place in the extreme conditions of the very early Universe and, 
indeed, they are the only clue to the origin of the Universe itself. Without 
observing and studying relic gravitational waves, there exists hardly any other 
way to test fundamental physical theories, such as superstrings and Af-theory. 
The understanding that the expansion rate of the very early Universe could 
be grossly different from the one governed by the radiatioDpdominated fluid 
exists for a long time, at least, from the study of Sakharov EH. In the recent 
years, however, there exists the unjust tendency to replace the notion and all 
the complexities of the very early Universe by a single word - "inflation" . The 
reminder of my contribution is devoted to relic gravitational waves and the 
associated issues of cosmological perturbations of quantum-mechanical origin. 

2 The inevitability of relic gravitational waves 

The theory of almost all astrophysical sources of gravitational waves is based 
on classical physics. In contrast, the theory of relic gravitational waves in- 
volves the elements of quantum physics. It is the power and beauty of relic 
gravitational waves that their existence relies only on general relativity and 
basic principles of quantum field theory, such as the concept of zero-point 
quantum oscillations. To understand the generation of relic gravitons it is 
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convenient to start from a classical picture. Imagine a classical gravitational 
wave propagating in a nonstationary universe. The nonlinear character of 
the Einstein equations provides the coupling of the gravitational wave field 
to the smooth gravitational field of the evolving universe. The gravitational 
wave itself can be arbitrarily weak, so that the quadratic terms of the wave 
field can be neglected. However, the nonlinearity of the Einstein equations 
leads to the presence of the time-dependent background gravitational field 
in the coefficients of the wave equation. The background gravitational field 
plays the role of the pump field-capable of increasing the amplitude of the 
propagating gravitational wave E3. In fact, the traveling wave is being ampli- 
fied, with the accompanying appearance of the wave traveling in the opposite 
direction. Together, they form an (almost) standing wave. This type of cou- 
pling of a wave-field to gravity is not something that is true for any wave-field. 
On the contrary, it is quite unusual. For instance, it does not take place for 
electromagnetic waves. 

The quantum element enters the picture when we realise that if the initial 
gravitational waves represent only the inevitable zero-point quantum oscilla- 
tions, the amplification mechanism should be at work. In precise technical 
terms, the quantum-mechanical Schrodinger evolution transforms the initial 
vacuum state of gravitational waives into a multi-quantum state known as 
squeezed vacuum quantum state It is the variance of phase that gets 
strongly diminished (squeezed) while the mean number of quanta and its 
variance get strongly increased. Of course, for this effect to be numerically 
considerable, the coupling of the wave field to the pump field should be con- 
siderable. The coupling is regulated by the ratio of the wave period to the 
characteristic time-scale of variations of the background gravitational field. If 
the pump field varies very slowly (adiabatically) , this ratio is very small, and 
practically nothing happens to the wave. In other words, during the adiabatic 
regime, the vacuum state remains the vacuum state. In the opposite (supera- 
diabatic) regime, the relative increase (as compared with adiabatic behaviour) 
of the wave amplitude, and the increase of the mean number of particles, is 
most dramatic. In the cosmological setting, the amplification mechanism is 
effective during the interval of time when the wavelength is comparable with, 
and longer than, the cosmological Hubble radius. Thus, the generating process 
(supepadiabatic (parametric) amplification of the zero-point quantum oscilla- 
tions E3) is a local and dynamical process, it does not rely on the existence 
of globally defined event horizons of the background space-time. Moreover, 
in cosmological models of practical interest, there is no horizons at all. It is 
also irrelevant which particular matter source drives the gravitational pump 
field, be it a scalar field or something else. In some of its basic aspects, the 
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quantum-mechanical generation of cosmological perturbations is similar to the 
laboratory-tested generation of squeezed light. 

The gravitational field of a homogeneous isotropic universe is given by 
the metric 

ds 2 = -c 2 dt 2 + a?(t)g ij dx i dx j = a 2 (r/) [-dry 2 + g ij dx i dx j ] , (1) 

where the scale factor a(t) (or a(r))) is driven by a matter distribution with 
some effective (in general, time-dependent) equation of state. The perturbed 
gravitational field can be written (for simplicity, in a spatially-flat universe) 
as 

ds 2 = a 2 {f 1 )[-dr 1 2 + (By + h ij )dx i dx j ], (2) 



(2tt) 3 / 2 



00 s=l,2 



(3) 



_±mx 



By expanding the functions hij(r], x) over spatial Fourier harmonics e 
we reduce the perturbed dynamical problem to the evolution of mode func- 

s s 

tions h n (v) f° r each mode n. Two polarisation tensors Pj_y(n), s = 1,2 have 
different forms depending on whether they represent gravitational waves, ro- 
tational perturbations, or density perturbations. In the case of gravitational 
waves, the polarisation tensors describe two polarisation states which are often 
called the "plus" and "cross" polarisations. 



For the quantized field, the quantities c n , c n 
operators satisfying the conditions 

st 



are annihilation and creation 



I c n, c m ] = <5 s / s c5 J (n - m) , c n |0) = 



(4) 



where |0) (for each n and s) is the fixed initial vacuum state defined at some 770 
in the very distant past, long before the superadiabatic regime for the given 

s 

mode has started. In that early era, the mode functions hn(v) behaved as 
oc e~ mT] , so that each mode n represented a strict traveling wave propagating 
in the direction of n. In the case of gravitational waves, the normalization 
constant C is V^Q^lpi- For cosmological density perturbations, which we will 
also discuss below, the normalisation constant is \p2Anlpi. 

The calculation of quantum-mechanical expectation values and correla- 
tion functions provides the link between quantum mechanics and macroscopic 
physics. Using the representation (|j|) and definitions above, one finds the 
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variance of metric perturbations: 



<0|M77,x)/^(r/,x)|0) = 




The quantity 

h?(n,r,) = ^n 2 ]T \L(v)\ 2 (6) 

s=l,2 

gives the mean-square value of the metric (gravitational field) perturbations 
in a logarithmic interval of n and is called the (dimensionless) power spec- 
trum. The power spectrum of metric perturbations is a quantity of great 
observational importance. It defines the temporal structure and amplitudes 
of the g.w. signal in the frequency bands of direct experimental searches. It is 
also crucial for calculations of anisotropy and polarisation induced in CMBR 
by relic gravitational waves and by other metric perturations. The inevitable 
squeezing makes the spectrum an oscillatory function of time, which can also 
be thought of as a consequence of the standing-wave pattern of the gener- 
ated field. At every fixed moment of time, the spectrum, as a function of the 
wave-number n, has many maxima and zeros. 

To find the power spectrum at any given moment of time (for instance, 
today or at the moment of decoupling of CMBR from the rest of matter) 
we need to know the mode functions at those moments of time. The mode 
functions are governed by the Hciscnberg equations of motion whose interac- 
tion Hamiltonian participates with the coupling function proportional to a' /a. 
These equations are, of course, equivalent to the perturbed Einstein equations. 
The second-order "master equation" describing gravitational waves is 113 

-0, (7) 

S 

where the functions are related to the mode functions hn(f]) by 

K(v) = a{v)L(v)- (8) 

One can view Eq. (^) as the equation for a parametrically disturbed oscil- 
lator (the term in brackets is the square of its variable frequency) , or as the 
Schrodinger equation for a particle moving in the presence of a potential bar- 
rier W(ri) = a" I a (while remembering that r\ is a time coordinate rather than 
a spatial coordinate). 
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As soon as the pump field (represented by 01(77)) is known, and since the 
initial conditions are fully determined, the mode functions are strictly calcu- 
lable. We know relatively well the behaviour of 0(77) at the matter-diminated 
(m) and radiation-dominated (e) stages of expansion. However, we do not 
know 01(77) at the most important precceding stage of evolution, which we 
call the initial (i) stage. Ideally, if the full evolution of 01(77) na d followed 
from some fundamental physical theory, the relic g.w. background would be 
uniquely determined, and we would have to live with what the fundamental 
theory dictates. In the absence of such a theory, we can only do calcula- 
tions for various possible models and compare the results with observations. 
Conversely, the observational restrictions on relic gravitons, or their direct 
detection, is our clue to the dynamics of the very early Universe. 

3 Calculating the power spectrum 

The scale factor 0(77) at the matter-dominated stage, governed by whatever 
matter with the effective equation of state p = 0, behaves as 0,(77) oc ?7 2 . It is 
convenient to write 01(77) in the explicit form 

a{i 1 ) = 2l H (i 1 -r lm f , (9) 

where Ijj is the Hubble radius today and rj m is a constant explained below. 
The moment of time "today" (in cosmological sense) is labeled by 77 = 77^ 
(the subscript R denoting "reception"). It is convenient to choose 

VR-Vm = I- (10) 

With this convention, a(rjji) = 21h, and the wave, of any physical nature, 
whose wavelength A today is equal to today's Hubble radius, carries the con- 
stant wavenumber rag = An. Longer waves have smaller n's and shorter waves 
have larger n's, according to the relationship n = 4itIh/X- For example, the 
ground-based gravitational wave detectors are most sensitive to frequencies 
around 30Hz — 3000Hz. The corresponding wavelengths have wavenumbers 
n somewhere in the interval 10 20 — 10 22 . 

The matter-dominated era was preceded by the radiation-dominated era 
with the scale factor 0(77) oc 77. Without any essential loss of generality, we 
assume that the transition from e era to m era was instantaneous and took 
place at some 77 = 772. The redshift of the transition is z eq : a{r}n) / '01(772) = 1 + 
z eq . It is believed that z eq is somewhere near 6 x 10 3 . In its turn, the radiation- 
dominated era was preceded by the initial era of expansion, whose nature and 
scale factor are, strictly speaking, unknown. To simplify the analysis, and 
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since the wave equations admit simple exact solutions in the case of power- 
law scale factors, we assume that the i era, similar to the e and m eras, was 
also described by a power-law scale factor. Following the early notations E3, 
we parametrize the i era by a(rf) oc |?7| 1+ ' 3 - The transition from i era to e era 
takes place at some 77 = rj\ and at redshift zf. a{rjn) / a{rj\) = 1 + z\. Further 
analysis shows (see below) that in order to get the right amplitude of the 
generated perturbations, the numerical value of Zi should be somewhere near 
10 30 . Whether the i era was governed by a scalar field ip (the central element 
of inflationary scenaria) or by something else, is irrelevant. What is relevant 
is the functional form of the gravitational pump field represented by the scale 
factor a{ji). 

We now write the full evolution of the growing scale factor explicitly: 

a(ri) = l \ri\ 1+f3 , V < Vi, '7i < 0, (3 < -1 , (11) 

a(r)) = l a e (r] ~ Ve)> m < V < V2 , (12) 

a(rj) = 2l H (r/ - r/ m ) 2 , m<V- ( 13 ) 

The continuous matching of 0(77) and a' (77) at the transition points determines 
all the participating constants in terms of Ih, Zi, z eq and (3. Thus, we are 
essentially left with the two unknown parametrs: Zi and (3. The Zi is primarily 
responsible for the amplitude of the generated perturbations, while the (3 is 
responsible for the spectral slope of the metric power spectrum. If (3 = —2, 
the generated spectrum is flat, that is, independent of the wave-number n. 
The flat spectrum is also called the Harrison-Zeldovich-Peebles spectrum and 
is characterised by the spectral index n = 1 . The relationship between the 
spectral index n and parameter (3 is n = 2/3 + 5. Inflationary models governed 
by scalar fields tp are uncapable of producing the "blue" power spectra, i.e. 
spectra with n > 1 ((3 > —2). Indeed, (3 > —2 requires the effective equation 
of state at the i stage to be e+p < 0, but this cannot be accomodated by any 
scalar field whatever the scalar field potential V(ip) may be 

For the CMBR calculations one also needs the redshift Zd ec of the last 
scattering surface r\ = t\e (with the subscript E denoting "emission"), where 
the CMBR photons have decoupled from rest of the matter: a(rjn) / a{rjE) = 
1 + Zdec- The numerical value of Zd ec is somewhere near 1000. 

Since the scale factor and initial conditions for the mode functions have 
been strictly defined, the power spectrum (^) is unambiguously calculable. 
We will present the results of calculation of h 2 (n,ri) taken at the time of 
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Figure 3. Plot of g 2 (x, 6 = 5) versus x. 



decoupling rj — t]e- The exact form of the power spectrum is 
h 2 (n,VE) = + z dec )n 4 \B\ 2 g 2 {x,b) , 



4< 



where 



(14) 



g 2 (x,b) = [p 1 {x)j 1 (bx) - p 2 (x)j- 2 (bx)] 2 , x = 



, b = 



_ 2 v /i+^ 



<■<] 



2y/l + Z eq y/T+Zde 



pi (x) = — [(Sx 2 — 1) sin x + Ax cos x + sin 3x] 



p2{%) = — r \(8x 2 ~ 1) cosx — 4x sin x + cos3a;l , 

and ji(bx), j-2(bx) are spherical Bessel functions. The quantity n A \B\ 2 is, in 
general, n-dependent and /3-dependent, but for the case of j3 = —2 it simplifies 
to 



nW = 4(1 + ^l P = -2. 



{l + z eq f 

In Fig. || we show the function g 2 (x, b) for 6 = 5. 



(15) 
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4 Detectability of relic gravitational waves 

We start from a discussion of manifestations aflrelic gravitational waves in 
the CMBR anisotropics. It can be shown (see tJ and references there) that 
the angular correlation function has the universal form 



ST ST 
^r0i)^7( e 2) 



0) =^^1(7^(008*), (16) 



4-7T 
1=2 



where Pi (cos S) are Legendre polynomials for the separation angle S between 
the unit vectors e! and e 2 , and the multipolc moments C; are explicitely 
determined by the metric mode functions. 

In Fig. ^ we show Eil by a solid line the graph of the function l(l + l)Ci cal- 
culated for the g.w. background with the power spectrum (f5). The cosmolog- 



ical parameters were taken, for illustration, as z eq = 10 4 , Zd ec — 10 3 , (3 = —2. 
The parameter zi is adjusted in such a manner [zi = 10 29 5 ) that the graph 
goes through the point 1(1 + l)Ci — 6.4 x 10~ 10 at I = 10, which agrees 
with observations. For comparison, the dashed line shows the same function, 
but calculated for the alternative (non-squeezed) background of gravitational 
waves of the same averaged power density. The remarkable (even if expected) 
result is that the background of non-squeezed (traveling) gravitational waves 
does not produce oscillations in the angular power spectrum Ci, whereas the 
background of squeezed (standing) gravitational waves does. The peaks and 
dips of the angular power spectrum are in close relationship with the maxima 
and zeros of the metric power spectrum of Fig. |^. It is unlikely that the peaks 
and dips caused by relic gravitational waves can be observationally revealed, 
but this calculation serves as a guidance for the related explanation (see be- 
low) of the actually observed peaks and dips i-\ that lay at a higher level of 
the signal. 

For purposes of direct experimental searches for relic gravitational waves, 
one needs to know the metric power spectrum in the present era. In Fig. [5] 
we show the piece-wise envelope of the r.m.s. quantity h(n, rjn) plotted as a 
function of frequency v in Hertz. 

This particular graph is derived under the assumption that the parameter 
(3 is /3 = —1.9, which corresponds to the spectral index n = 1.2. A few years 
ago (see, for example, t£l), this value of n was quoted as the one that follows 
from the lower-order multipoles of the CMBR data. The most recent analyses, 
which include data on the large scale structure, seem to prefer a smaller value 
of n. Although the often cited central values of n are still somewhat larger 
than 1, the la error-bars do not usually include the value n = 1.2. One should 
remember, however, that the current analyses of the data are consistent only if 
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Figure 4. The solid line depicts the plot of 1(1 + 1)C; versus I, normalized such that at 
I = 10, we have 1(1 + 1)C; = 6.4 X 10 — 10 , which tallies with observations. The dashed line 
is the corresponding plot for the non-squeezed g.w. background. We take fi = —2, and the 
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10 4 and z^e 



10 3 . 



one assumes that about 70 percent of "matter" in the present Universe resides 
in the cosmological A term or in a "dark energy/quintessence" , which is not 
subject to gravitational clustering. It is still not clear whether this conclusion 
is one of the greatest discoveries about the Universe, or it may be a result of 
improper evaluation of the role of metric perturbations in the CMBR analysis. 
We will comment on this issue in the next section. In addition, one should 
bear in mind that the flat (n = 1) and "red" (n < 1) spectra possess a 
theoretical difficulty: the integral over wavenumbers n for the mean square 
value of the metric (gravitational field) perturbations diverges at the lower 
limit n — ► 0, i.e. for very long waves. The removal of these divergencies would 
require some extra assumptions. In any case, we will discuss the prospects of 
direct detection of relic gravitational waves assuming that the graph in Fig. 
^ is correct or close to correct. 

According to analytical formulas and Fig. ^, the theoretically expected 
r.m.s. value at v = 100Hz is h th = 10 -25 . The target value of the LIGO-IIB 
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Figure 5. Expected spectrum h(v) for the case /3 = —1.9. 



at this frequency is h ex — 10~ 23 . The remaining gap of 2 orders of magnitude 
can be covered by the cross-correlation of sufficiently long strips of data from 
two or more detectors. The S/N ratio will be better than 1, if the common 
integration time exceeds r = 10 6 sec. This does not look like a hopeless task. 
The situation is somewhat better in the LISA frequency interval. According 
to evaluations of Ref. a 1-year observation time can provide a S/N ratio 
around 3 in the frequency interval (2 x 10~ 3 — 10~ 2 )Hz. These estimates do 
not include a possible improvement in the S/N, which can be achieved if one 
succeeds in finding a sophisticated method of exploiting, at these frequencies, 
the ever-present squeezing. This improvement would certainly be possible in 
a very narrow-band detector, but it is unclear whether it can ever_be realised 
in a realistic broad-band detector (for a discussion, see E3 and Q). At the 
present time, it is difficult to say at which scales the relic gravitational waves 
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will be first discovered. This can happen in direct experiments at relatively 
small scales, or, indirectly, at cosmological scales, through the anisotropy and 
polarisation measurements of the CMBR. 

5 Relic gravitational waves and primordial density 
perturbations 

The superadiabatic (parametric) mechanism is always operative for gravita- 
tional waves, but its applicabilty to density perturbations requires special as- 
sumptions about properties of matter. Of course, the notion of density pertur- 
bations includes the associated metric perturbations of Eq. |^; the matter 
and gravitational field (metric) variables are always linked by the perturbed 
Einstein equations. If the i stage was governed by a scalar field ip coupled 
to gravity in a special manner (a hypothesis usually assumed in inflationary 
scenaria), the quantum-mechanical generation of primordial density pertur- 
bations, in addition to the inevitable generation of relic gravitational waves, 
becomes possible. The underlying physical theory for these two respective 
fields is exactly the same, and therefore the general properties of relic gravita- 
tional waves and primordial density perturbations (squeezing, standing-wave 
pattern, spectral slopes, numerical values of the amplitudes, etc.) are al- 
most identical. In particular, relic gravitational waves and primordial density 
perturbations should provide approximately equal contributions to the lower- 
order multiplies (starting with the quadrupole moment Cq) of the CMBR 
anisotropies E3. Some differences arise, however, at the late stage of their 
evolution. In the matter-dominated era, and for wavelengths comfortably 
shorter than the Hubble radius, numerical values of the metric components 
associated with density perturbations grow in proportion to a(r\) (gravita- 
tional instability), whereas the amplitudes of gravitational waves decrease in 
proportion to 1/0(77) (adiabatic behaviour). This explains why the metric 
amplitudes of these two fields, being of the same order of magnitude at scales 
comparable with the prersent-day Hubble radius, are significantly different at 
much smaller scales. This difference makes gravitational waves subdominant 
in their contribution to the CMBR multipoles near the peak at I ~ 200, and 
makes them so small and difficult to detect at laboratory scales. 

As a benchmark for the spectral index one can take (3 — —2, i.e. flat 
spectrum. The CMBR anisotropies caused by gravitational waves are shown 
in Fig. The lowest excited multipole is I = 2 (quadrupole). In contrast, 
density perturbations produce all multipoles, including I = (monopole) and 
I = 1 (dipole). In general, CMBR anisotropies caused by density pertur- 
bations have contributions from 3 different ingredients (for a recent detailed 
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discussion, seell). Specifically, the ingredients are: the accompanying metric 
perturbations (mostly driven by the gravitationally dominant matter, which 
is, according to the existing views, a presureless cold dark matter (CDM)), 
the intrinsic temperature variations at the last scattering surface, and the ve- 
locities (with respect to the coordianate system comoving with the CDM) of 
the last scattering electrons. 

The angular power spectrum l(l + l)Ci is usually presented as an outcome 
of numerical calculations, where all contributions are being mixed up. How- 
ever, the evolutionary equations, on which the numerical codes are based, 
allow analytical calculations for the lower multipoles, i.e. for l's up to 
I w (30 — 40). It is often claimed that the analytical limit of the employed 
evolutionary equations produces the lower-order multipoles C; according to 
the formula 

This formula predicts an infinitely large monopole moment Co, and the dipole 
moment C\ only a factor 3 greater than the quadrupole moment C2; both 
predictions are in a severe conflict with observations. Moreover, this formula 
misses the quite steep growth, beginning from I = 2, of the function 1(1 + l)Ci, 
which must take place for purely gravitational reasons, that is, when the 
metric perturbations are treated correctly EJ. 

It was demonstrated til that formula [l7] is incorrect and cannot be the 
outcome of correct equations. Specifically, formula ^ arises only if one unjus- 
tifiably—neglects 3 out of 4 terms in the exact formula (43) of the Sachs- Wolfe 
paper E3. The correct handling of metric perturbations makes the monopole 
moment Co finite and small, whereas the dipole moment C\ becomes 5 orders 
of magnitude greater than the number following from formula and, thus, 
the theoretical (statistical) C\ comes in agreement with the observed value of 
C\. It was also shown that the growth of the function 1(1 + 1)C would be 
unlimited, if one were allowed to extrapolate the flat spectrum to arbitrarily 
large wave-numbers n. This is illustrated in Fig ^|. 

In reality, however, the modulating (transfer) function of the metric per- 
turbations bends the angular spectrum down and introduces a series of peaks 
and dips, analogous to the peaks and dips caused by gravitational waves in 
Fig. |j. We will discuss peaks and dips shortly, but first we shall say a few 
words about the importance of theoretical calculations of the expected (sta- 
tistical) value of the dipole moment C%. Because of statistical properties of 
the generated perturbations, all statistical multipoles C are independent of 
the observer's position. 
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Figure £L_The multipole distributions as described by the exact formula (solid line) and 
formula hi (dotted line) for the flat primordial spectrum (3 = —2 (n = 1). 



The derivation of formula [17| takes into account the perturbations of all 
wavelengths. The integral over n, which gives rise to this formula, is formally 
extended from n = to arbitrarily large n's. If this formula were correct, it 
would be applicable all over the linear regime of cosmological perturbations 
and could possibly break down only due to nonlinearities at relatively short 
scales. The gap in 5 orders of magnitude between the statistical prediction 
of this formula and the actually observed C\ would need to be attributed 
either to extremely improbable realisation of the perturbed gravitational field, 
or to the quite local conditions where nonlinearities set it and where the 
linear approximation is not valid. In the latter case, numerical values of Ci's 
observed by inhabitants of the not-too-distant galaxies would not need to 
have anything in common with-each other and with the C\ observed by us. In 
contrast, calculations of Ref. Ell show that numerical value of the statistical 
C\ is most sensitive to scales where the metric power spectrum bends down. 
In other words, the corresponding integral is primarily accumulated at scales 
of the order of (100 — 2Q0)Mpc, i.e. well within the linear regime. This result 
comes about from the "gradient of gravitational potential" terms (see also ta) , 
which are totally neglected in the derivation of formula This means that 
the galaxies within, at least, the radius of (100 — 2QQ)Mpc, and, at least, the 
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field galaxies, are expected to have approximately equal observed Ci's. This 
number should also be of approximately the same value as our own C\ (as soon 
as our measured C± turns out to be close to the statistical C\). In principle, 
this picture is observationally testable. Of course, we cannot place ourselves 
at different galaxies and measure the C\ there. But, it is not excluded that 
the CMBR dipoles seen at other galaxies can be evaluated indirectly, possibly 
with the help of the Sunyaev-Zeldovich effect. It would be very interesting 
if these observational evaluations could be done. In any case, the calculation 
of the theoretical (statistical) C\ is as important as the calculation of other 
multipoles. It would be embarassing if some of currently popular cosmological 
models would predict the statistical C\ way out of what is actually observed 
by us and by observers in other galaxies (assuming that we can evaluate their 
dipoles). 

It is argued in Ref. EJ that the (squeezed) metric perturbations associated 
with density perturbations are mostly responsible for the recently observed 
peaks and dips in the angular power spectrum 1(1 + 1)C;. This interpretation 
is different from the concept of "acoustic peaks" based on the existence of 
the plasma sound waves at the last scattering surface. The rigorous evolution 
of the quantum-mechanically generated density perturbations through the 
radiation-dominated and matter-dominated eras, allows one to find the metric 
power spectrum, defined by Eq. ^, at the last scattering surface r\ = The 
analytical formula for h 2 (n,r]E) is given by 



C 2 n 4 \B\ 2 (l + z eq ) (smy 2 \ 2 (300 - 20p 2 y 2 + p\|) 
2ir 2 A8l 2 H \ V2 ) 200 



h{n,7i E ) = ^ — 2 [ — ) ^ , (18) 

where 



2/2 = - 7= n =, and p = 



2V3-v/l 



2V3^/1 + Z eq ' VI + z dec 

The graph of the oscillating function 

f(x,p)= 2 [300 - 20p 2 x 2 + p 4 x 4 ] , 

where x = t/2 and, for illustration, we take p — 8, is shown in Fig. ^. 

The oscillations in the metric power spectrum are a direct consequence 
of squeezing (standing-wave pattern) of the primordial density perturbations. 
These oscillations find their reflection in the peaks and dips of the angular 
power spectrum. The main difference with the concept of "acoustic peaks" 
lays in the two facts: the amplitudes of metric perturbations are larger than 
the amplitudes of matter density variations in the photon-electron-baryon 
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Figure 7. The plot depicted by the solid line is that of f 2 (x,p = 8) versus x. The dashed 
line shows the behavior of the function (sinx/x) 2 X 10 6 . 



fluid, and the wave-number periodicity of metric perturbations is governed 
by a sine frrnction, instead of a cosine function for acoustic perurbations. It 
is argued Eil that the actually observed features are better described by the 
derived analytical formulas than by the concept of "acoustic peaks" . 

The study of density perturbations is directly relevant to the problem of 
relic gravitational waves. If cosmological perturbations are indeed generated 
quantum- mechanically (what the present author believes is likely to be true) , 
many properties of primordial density perturbations and relic gravitational 
waves should be common. By studying density perturbations (for instance, 
by proving the presence of squeezing) we will be more confident about the 
properties of relic gravitational waves. 

One should be aware that the story of quantum-mechanically generated 
cosmological perturbations is dramatically different in inflationary scenario. 
Inflationists love to speak about "stretching the quantum fluctuations of the 
inflaton to cosmological scales" . As a result of this "stretching" , they manage 
to produce the arbitrarily large amplitudes of today's density perturbations. 
Indeed, the "standard result" of inflationary scenario predicts the infinitely 
large density perturbations, in the limit of the flat spectrum (n = 1), through 
the set of evaluations: Sp/p ~ h S ~ C ~ H 2 /(p ~ V 3 / 2 (<p)/V' {(p) ~ - n. 

[The "quantization", as understood by inflationists, of the necessarily cou- 
pled system of scalar field and metric perturbations is such, that the so-called 
Bardeen's gauge-invariant quantity £ is arbitrarily large already at the begin- 
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ning of the amplifying (superadiabatic) regime, and then, being "conserved" , 
this arbitrarily large number translates to the end of the amplifying regime, 
and to the present time.] Of course, this "standard result" is in full disagree- 
ment not only with the theoretical quantum mechanics, but with available 
observations too, as long as the error-boxes of the observationally derived 
spectral index n are centered at n w 1 and include the value n = 1. How- 
ever, by composing the ratio of the gravitational wave amplitude hx to the 
predicted divergent amplitude of the scalar metric perturbations h$ (the so 
called "consistency relation": hx/hs ~ yl — n), inflationary theorists substi- 
tute their prediction of arbitrarily large density perturbations for the claim 
that it is the amount of gravitational waves that should be zero, or almost 
zero, at cosmological scales and, hence, down to laboratory scales. This claim 
has led to many years of mistreatment of a possible g.w. contribution to the 
CMBR data. The "standard inflationary result" is maintained by inflationists 
and their followers until now. For instance, the recent review paper claims 
that the initial spectrum of gravitational waves is "constrained to be small 
compared with_thc initial density spectrum" . ft is only in a few recent papers 
(for example, Ej) that the inflationary "consistency relation" is not being used 
when analyzing the CMBR and large scale structure observations, with some 
interesting conclusions. The inflationary claim of the "negligible" amount 
of gravitational waves influences also the LISA community. It is sometimes 
opined that the LISA parameters should be better optimised for detection 
of g.w. noise from the multitude of binary white dwarfs, rather than for a 
search for relic gravitational waves. There seems to be some logic behind this 
opinion. If one hears so often about "excellent agreement" of inflationary 
predictions with observations, and if the central of these predictions is "neg- 
ligible" amount of relic gravitational waves, there is no need to bother about 
them. Many of the contemporary writers manage to do both: to praise infla- 
tionary predictions and to speculate about "testing inflation" with the help 
of gravitational waves, including the CMBR polarisation measurements and 
post-LISA missions. Apparently, they plan to do this at the signal to noise 
ratio level of 10 -20 , or so, for "inflation predicted" gravitational waves. For 
some reason, those authors do not suggest to test inflationary predictions right 
from the predicted divergent amplitudes of density perturbations. Sometimes, 
inflationary calculations based on the "standard result" and "consistency re- 
lation" arrive at gravitational-wave amplitudes which are only 1-2 orders of 
magnitude smaller than, or even comparable (for the spectral index n suffi- 
ciently smaller than 1) with the metric amplitudes of density perturbations. 
Of course, this does not make the fundamentally wrong "standard result" 
a correct one. One can find more reading about the "standard inflationary 
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result" and its critique, including the earlier claim of inflationists that their 
divergent amplitudes of density perturbations are justified by the (incorrect 
and now abandoned) .concept of "big amplification during reheating", in the 
end of Sec. 6 of Ref. □ and in references there. 



6 Conclusions 



It is important to understand that the existence of relic gavitational waves 
relies only on general relativity and basic principles of quantum field theory. 
The detection of relic gravitational waves will be a direct probe of the expan- 
sion rate of the very early Universe, and not a test of some extra assumptions, 
such as the dominance of a scalar field or the form of its potential V ((f). If the 
observed large-scale CMBR anisotropy ST/T is indeed caused by cosmological 
perturbations of quantum-mechanical origin, the contribution of relic gravi- 
tational waves to 5T/T must be considerable. The extrapolation of the g.w. 
metric power spectrum to shorter wavelengths shows that, depending on the 
still existent uncertainty in the spectral index, relic gravitational waves may 
be measurable by the advanced ground-based or space-based laser interferom- 
eters. The ever-present squeezing must manifest itself in the periodic structure 
of the metric power spectra and, at cosmological scales, in the oscillatory be- 
haviour of the CMBR multipoles Ci as a function of I. When evaluating the 
scientific importance and prospects of discovery of relic gravitational waves, 
it is necessary to respect the conclusions of quantum mechanics and general 
relativity, and not the statements of popular fallacies. 
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